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I. INTRODUCTION 

In the mid 1980's Harris, Lubensky and coworkers 
(HL) [l|, 13, "I developed a seminal field theoretic model 
for the unified description of percolating random resis- 
tor networks (RRNs) and a class of diluted continuous 
spin systems (such as the rr-y-model). Their approach, 
based on ideas by Stephen Q, turned out to be very 
fruitful. In the course of the years, it provided a foun- 
dation for the exploration of various criticalproperties, 
not only of the s-y-model and the RRN J3,J3, la Q, but 
also of random resistor diode networks 0, Oj , the swiss 
cheese model llfil. fTT| and random networks of Josephson 
junctions [l2L Il3|. Moreover, it fostered the computa- 
tion of several fractal dimensions of isotropic 0, 0, 
and directed 0, 0] percolation clusters and aided study- 
ing multifractality in isotropic [Isl ITfll Eol ] and di- 
rected |2lJ, |22| percolation. Also, the HL model helped 
to improve the understanding of the vulcanization tran- 
sition 23]. 

An important role in the HL model is played by the 
exponent 0. It describes the power law behavior of the 
average resistance Mjj between two connected points x 
and x' at the percolation threshold, 



M 



(i) 



r'\<t>/v 



(1.1) 



Here, v is the critical exponent of the correlation length 
£ ~ \P — Vc\~ u where p is the probability that controls 
the dilution of the network and p — p c marks the critical 
point. 

Originally [2j <f> was thought to be the first member, 
<j) = <f>i, of an entire family {<j) m , m = 1,2,...} of expo- 
nents with the mth member describing the mth cumulant 
of the resistance, 

M { ™ ] ~ \x-x'\^ /u . (1.2) 

The (j) m were calculated to one- loop order Q. Shortly 
later Rammal, Lemieux, and Tremblay |f| showed that 



the higher crossover exponents 
that higher cumulants scale as 



,1 are irrelevant and 
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nm<j>/u 
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i.e., they display so called gap scaling. Furthermore, 
Rammal, Lemieux, and Tremblay argued that the 4> m >i 
are important for corrections to scaling and specifically 
that these are governed by exponents u> m = {m<j)~(j) m )/v. 
For example, the average resistance was thought to be- 
have as 



m 

and the conductivity £ as 



\x — x 
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(1.4) 



(1.5) 



and 



Here, t is the conductivity exponent t = (d— 2)v- 
A m and B m are non- universal amplitudes. 

In this paper we take up the issue of corrections to 
scaling in RRN anew. Our careful analysis reveals that 
previous work on this subject is erroneous and that the 
4>m>i are meaningless, at least as far as corrections to 
scaling in RRN are concerned. It turns out that the pre- 
vious studies overlooked a crucial feature of irrelevant 
field theoretic operators, viz. that they tend to mix un- 
der renormalization. Taking this subtlety into account, 
we calculate corrections to scaling for the average resis- 
tance and the conductivity. Moreover, we re-examine 
corrections to scaling in continuous spin models. Due to 
a difference in the symmetry properties that plays no role 
for the leading scaling behavior, the corrections to scaling 
in continuous spin models and RRN are somewhat dif- 
ferent. We determine the correction to scaling exponents 
for the first cumulant of the spin orientations. 

The outline of our paper is as follows. In Sec. ITT1 we 
present a few basics about the RRN and the s-y-model. 
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Then we state the field theoretic Hamiltonian defining 
the HL model and sketch its physical contents. Next we 
collect the irrelevant field theoretic operators that lead to 
the corrections to scaling we are interested in. These op- 
erators can be classified into 2 groups, viz. specific opera- 
tors and general percolation operators. Section lHfl brieflv 
reviews known field theoretic results for the HL model 
including the leading scaling behavior and the leading 
correction to it. In Sec. II VI we analyze the equation of 
motion implicit in the HL Hamiltonian. This analysis 
provides us with several results that are valid to arbi- 
trary order in perturbation theory. Section Ivl comprises 
our renormalization group analysis of the specific oper- 
ators. The general percolation operators are scrutinized 
m Sec. El In Sec. EH we derive our final results for the 
critical behavior of the average resistance etc. The main 
part of our paper concludes with several remarks given 
in Sec. IVIIll There are 3 Appendices. In Appendix lAl 
we derive two useful identities that help us exploiting 
the equation of motion. Appendix iBl features some gen- 
eral considerations on composite fields in the HL model. 
Appendix|Ocontains details on the computation of Feyn- 
man diagrams. 



II. THE HARRIS-LUBENSKY MODEL 

A. Random resistor networks and the diluted 
x-y-model 

Both the RRN and the diluted x-y-model can be de- 
scribed by a Hamiltonian of the type 

H = -J2 U *A^-^) , (2-1) 

(i,3) 

where t?j is a continuous dynamical variable pertaining 
to a site i of a rf-dimensional hypercubic lattice and the 
sum runs over all nearest neighbor pairs on this lattice. 
For the RRN, corresponds to the voltage Vi at site i 
and is defined on the interval [—00, 00]. 

Ui, j {V) = -\a i , j V 2 (2.2) 

is the electrical power dissipated on the bond between i 
and j, with aij denoting the conductance of this bond. 
(Tjj is a random variable that takes on the values a and 
with probability p and 1 — p, respectively. In the case 
of the x-y-model, i?j is the angle ifi that specifies the 
orientation of the spin at site i and is defined on the 
interval [— 7T, tt]. Here, 

U it j (ip) = K itj cos ip (2.3) 

with Ki j being the exchange integral. In the diluted x- 
y-model the exchange integral is assumed to take on the 
values K and respectively with probability p and 1 — p. 



B. The Harris-Lubensky Hamiltonian 

Based on the Hamiltonian 12.1fl HL derived a field the- 
oretic model that can be written as 

h = /^{E(r 2 +^ Vs ) 2 +f( v ^) 2 +r 3 ) 



i J 

For details on the derivation we refer the reader to 
Ref. @]. The order parameter field s(x, 9) lives on a con- 
tinuous d-dimensional space with the coordinates x. It 
is subject to the constraint 

X>M) = 0. (2.5) 
e 

The variable 9 is a replicated analog of the original dy- 
namical variable d and takes on discrete values on a D- 
dimensional torus (the replica space) . Formally, 9 is given 
by 9 — kA9, where A: is a D-dimensional vector with in- 
teger components k^ and —M < < M. To recover 
the physical situation, on has to take the replica limit 
D -> 0, M -► 00 with M D -► 1 and A9 = Oq/VM. 
In this limit, 9q plays the role of a redundant scaling 
parameter, i.e., the theory is independent of its value. 
The parameter r is proportional to p c — p, i.e., it spec- 
ifies the distance from the critical point, w is propor- 
tional to cr _1 or K^ 1 , respectively. The A4, finally, are 
irrelevant field theoretic operators (rotationally invariant 
monomials constructed from the fundamental field s and 
its derivatives in real and replica space). Specifics of the 
irrelevant operators will be given below. 

For w = fi = 0, the HL Hamiltonian describes the N- 
state Potts model with N = J2g 1= ( 2M ) D ■ In this case 
we have Sjv, the group of all permutations of N objects, 
as the internal symmetry group. If w ^ 0, this symmetry 
is reduced to O(D), the group of orthogonal rotations in 
the replica space. A particular scaling symmetry of the 
Hamiltonian will be important as we go along. Namely, 
H is invariant under the rescaling w — > bw because the 
scaling factor b can be absorbed into the redundant pa- 
rameter 9q. 

In the following we use that J2g ■ • ■ ~ (A6*) D J2s ■ • • ~ 
J d D 9 . . . and abbreviate the latter integral by Jg---- The 
approximations involved here become exact in the replica 
limit. 

C. Physical contents 

To fully appreciate the physical contents of the HL 
model it is helpful to consider the replica space Fourier 
transform 

= J exp 0) s(x,e) (2.6) 
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of the order parameter. For completeness we mention 
here that the constraint i|2.5[l translates upon replica 
space Fourier transform into 



lfo(x)=0. 



(2.7) 



This constraint is intuitively clear because i/'g ( x ) is a con- 
stant and hence does no qualify as an order parameter. It 
will play an important role in our renormalization group 
(RG) analysis. 

The value of the quantity V'a( x ) * s that its correlation 
functions 



G(x,x';A 



(2- 



provide convenient access to physical observables. In the 
case of the RRN one has 



G 



(x,x';X) 



= ( exp 



-yfl(x,x') 



(2.9) 



where i?(x,x') = ([V x — V x >] 2 ) is the macroscopic resis- 
tance between the points x and x' and (• • • )c denotes 
averaging over all configurations of the diluted lattice. 
Hence, G(x,x';A) is a generating function for the mo- 
ments of the resistance distribution 

A4 i) (x,x') = <i?(x,x') Z )' c , (2-10) 

where the prime indicates averaging subject to the con- 
straint that x and x' are connected. 

For the diluted x-y-model the Hamiltonian (|2.1|) is not 
Gaussian and one is led to 



(2.11) 




-I ' G 



where (• • -p c ' stands for the cumulants with respect to 
the average (•••) and Ki{\) = ^f =1 [A(°>] 21 . Thus, 
G(x, x'; A) is a generating function for the cumulants 



G«(x,x') = 



_fx - <fx 



(2.12) 



which measure the fluctuation of the angular variables. 
c[p \ in particular, is related to the spin- wave stiffness. 



D. The irrelevant operators Ai 

The main goal of this paper is to analyze corrections 
to scaling. The leading correction to scaling for the aver- 
age resistance etc. is well known since it is described by 
the so-called Wegner exponent, see Sec. MI Bl below. The 
next to leading corrections stem from irrelevant operators 
which scale as Ai ~ n s at the upper critical dimension 



d c = 6, where \x is some inverse length scale. These op- 
erators can be classified into two groups. The first group 
consist of operators having at least two derivatives with 
respect to the replica space. The operators belonging to 
this group will be referred to as specific operators. In the 
case of the RRN, the specific operators are given by the 
0(D) invariant composite fields 



Ai 
A 2 
A 3 



f v 2 ( / ,v r ;, 



sV 2 V?,s. 



G 



(2.13a) 
(2.13b) 
(2.13c) 
(2.13d) 



Each replica space derivative is accompanied by a factor 
^/w to ensure invariance under the rescaling w — > bw. 

When considering the x-y-model, we have to admit all 
of the hypercubic invariants of Vg rather than of V 2 ,. 
Thus, we are confronted with a further operator with the 
naive dimension 8, namely 



•Ac 



v 2 ( 



(2.14) 



Note that this operator breaks the 0(£>)-symmetry in 
replica space. 

As we go along we will see that the family 
{^o ; • • • j ^3} is associated with corrections to scaling of 
the average resistance and conductivity in RRNs. Hence, 
we refer to this family as resistor specific. The larger 
family {A c , Ao, ■ • ■ , ^3} is associated with corrections of 
the spin orientation cumulants in diluted continuous spin 
systems. Thus, we say this family is spin specific. 

The next group of operators leads to corrections for 
both pure percolation and the respective specific behav- 
ior. We refer to these operators as general percolation 
operators. The pure percolation behavior results from 
the HL model in the limit w = 0. Thus, the second 
group involves only operators without derivatives with 
respect to the replica space. They are given by the Sn 
invariant composite fields 



Ai 

A:-, 

A> 
A 7 



'V 2 s 



(2.15a) 
(2.15b) 
(2.15c) 
(2.15d) 



At first sight, three additional general percolation opera- 
tors with naive dimension eight seem to matter, namely 
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Ag = 


2 V J 


v-72 

SV S , 

8 


(2.16b) 


Aio = 




Je 


(2.16c) 


An = 




8 


(2.16d) 



However, these operators can be neglected in calculat- 
ing corrections to scaling. Upon Fourier transformation 
of A4 and As one sees that these operators coincide for 
vanishing external momentum. Thus, it is sufficient foe 
our purposes to keep A4. Ag to An are merely total 
derivatives of lower dimensional operators. Their Fourier 
transformed counterparts are proportional to some (pos- 
itive and even) power of a external momentum and hence 
these operators cannot contribute to a translational in- 
variant Hamiltonian. This situation is similar for Ag ■ We 
choose, though, to keep Aq in our analysis for 2 reasons: 
(i) To exemplify explicitly that this kind of operator does 
in the end not contribute to the corrections to scaling, (ii) 
Based on the scaling symmetries of TC we can draw exact 
conclusions on Aq which then can be compared to the 
results of our explicit 1-loop calculation, i.e., retaining 
Aq allows for a later consistency check. 



III. A BRIEF REVIEW OF KNOWN 
RENORMALIZATION GROUP RESULTS 



G e = T(l + e/2)/{4n) d / 2 is a dimension-dependent nu- 
merical factor. We have determined the renormalization 
constants Z...(u) to second order in u in a 2-loop calcula- 
tion. From the scheme (|3.1[) we derived a Gell-Mann-Low 
RG equation for the vertex functions with n ampu- 
tated external legs, 



d a d 



d . d 

TK— h WC, - 



dfj, du dr dw 
x T (n) (\x,x\]T,u,w,ti) =0, 



(3.2) 



with the Wilson functions f3 (u) = 
C(u) = M ^|o and 7 (u) 



, , d In r , 



0§*|o, «(«) 



Solving Eq. (|3.2I) at the infrared stable fixed point u* 
leads to the scaling behavior of the vertex functions [2(| , 

({x, a} ; r, U , w, M ) = r ( d - 2 +")"/ 2 
x (Ux, a} ; r 1 ' v t, u* , r*/ v w, m) ■ (3.3) 

The exponents rj = 7*, v — (2 — k*), and <p = v{2 — 
£*) [7* = 7(m*), k* = k(u*) and so on] arc the usual 
critical exponents for percolation and the RRN. In the 
present paper we work to 1-loop order. To this order, 
u* = 2e/7 + 0(e 2 ) and the critical exponents are given 
by ry = -e/21 + (9(e 2 ), v = 1/2 + 5e/84 + 0(e 2 ) and 
4>= l + e/42 + (9(e 2 ). 

From (|3.3|) one can extract the leading scaling behavior 
of various observables. Exploiting Eq. (|2.9() . for example, 
it is straightforward to derive Eq. (|1.3|) for . 



Many of the critical properties of the HL model are 
well known. Substantial contributions stem, inter alia, 
from HL. In earlier work we have investigated the HL 
by using the powerful methods of renormalized field the- 
ory. Here we briefly review parts of this work to provide 
background and to establish notation. 

A. Renormalization and Scaling 

In Ref. |6( we have studied the HL model with all the 
fi equal to zero. In particular we have determined the 
renormalizations 

s^s = Z 1/2 s, (3.1a) 

r^f = Z" 1 Z T T, (3.1b) 

w — > w = Z~ 1 Z w w , (3.1c) 

g->g = Z- 3 / 2 Z g9l (3.1d) 

u = G e [T e g 2 (3.1e) 

in dimensional regularization and minimal subtraction 
of e-poles [24|. Here, the denotes bare, unrenormalized 
quantities. The coupling constant u is introduced be- 
cause it is convenient and dimensionless. e — 6 — d mea- 
sures the deviation from the upper critical dimension. 



B. The leading correction to scaling 

The leading correction to scaling is, as usual, governed 
by the so-called Wegner exponent. This leading correc- 
tion emerges when the renormalized coupling u is not 
exactly equal u* since the renormalization flow has not 
arrived at its fixed point yet. Such a case occurs typically 
when there is a finite momentum cutoff reminiscent of a 
non- vanishing lattice spacing. 

Taking the leading correction into account, the scaling 
behavior, e.g., of the average resistance (|1.1|) becomes 

m£> ~ \x - x'\*l v [1 + A\x - x']-"] , (3.4) 

where A is a non- universal amplitude and lu = (3' (u*) is 
the Wegner exponent, lu can be calculated without much 
effort to third order in e upon using the three-loop result 
for (3 (u) obtained by de Alcantara Bonfim et al. [23 • To 
the order we are working here, the Wegner exponent is 
given by 

uj = e + 0(e 2 ) . (3.5) 

Expression similar to (|3.4|) hold for the conductivity and 
so on. 
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IV. CONSEQUENCES OF THE EQUATION OF 
MOTION 



Several of the correction to scaling exponents origi- 
nating from the specific and the general percolation op- 
erators can be derived without resorting to an explicit 
perturbation calculation. We will do so by analyzing the 
classical equation of motion. The so obtained results have 
the virtue of being rigorous in the sense that they hold 
to arbitrary in perturbation theory. In the remainder we 
will frequently encounter the so-called scaling dimension 
x^ of an eigenoperator A. Formally, x^ is defined via 
the rescaling -4.(x) = £ XA A(£yi) when A is inserted into 
vertex or correlation functions at the critical point. 

Before we turn to the consequences of the equation of 
motion we take a little detour and extract a consequence 
of the invariance of TL under the rescaling w — > bw. It 
follows from this rescaling invariance that, in the case 
w 7^ 0, the operator w /g(Vgs) 2 is marginal and its scal- 
ing dimension is equal to d. Note that A$ results from 
the application of the Laplacian V 2 to the this marginal 
operator. Consequently the scaling dimension of Aq is 
given by 



x Ao 



d+2. 



(4.1) 



Now we turn to the equation of motion. It is well 
known that for every independent equation, that fol- 
lows from the equation of motion, there is an eigenop- 
erator of the RG with a scaling dimension that can be 
expressed in terms of the scaling dimensions of lower- 
dimensional operators. For some background on four- 
and six-dimensional operators in the HL model we refer 
to Appendix iBl 

The classical equation of motion derived from the 
Hamiltonian l|2.4|l and the constraint i|2.5[l reads 



_ SH 

Ss f=o 
-V 2 s - wVls 



2 V N 



0. (4.2) 



In the following we work in the limit D — > 0, i.e., 
N = (2AI) D — > 1. We consider the lower dimensional 
operators 



Bi = w\7 2 g s 



B 2 = V 2 



B 3 =H'. 



(4.3) 



The scaling dimensions of the first two follow from their 
renormalizations as 



d-2 + r\ <j) 



XBi = ^ h — , x B2 = 

The scaling dimension of £> 3 , 

x Bi = d-2 + rj B3 



d + 2 + rj 



(4.4) 



(4.5) 



is not the same as the scaling dimension of the operator 
Bo = J g s 2 . Bq belongs to the trivial representation of the 



permutation symmetry group Sn, whereas the operators 
Bi with i = 1, 2,3 transform like s(x, 9), i.e., they belong 
to the fundamental representation of Sn [permutations 
of N objects 4>i with constraint ^ i=1 <fii — 0]- Hence, 
r)B 3 remains to be determined. In Appendix IaI we derive 
a Ward identity, Eq. I)A17(I , that implies 



VB 3 = 



(4.6) 



to arbitrary order in e expansion. Consequently, we ob- 
tain the exact result 



d + 2 



x Ba 



n 



(4.7) 



Now we can extract the actual consequences of the 
equation of motion for our prime targets of interest, viz. 
the eight-dimensional operators. In Appendix IaI we show 
that 



»(xi, 00 • • • s( Xn ,9 n ) T ■ W'(x)) = ^<5(x - X,) 



x ^s(xxA) • ■•s(x i _ 1 ,0 i _i)^ : '(x i) 0i) 
x s(xi+i, 9 i+1 ) ■ ■ ■ s(x„, 4)) ■ (4.8) 
where T is some composite field and where we defined 
■F-H'(x) = [ F(x,9)H'(x,9). (4.9) 



Upon renormalization, on finds that the identity l|4.8() 
implies the scaling relation 



d - X s + Xjr 



d + 2- 



V 



(4.10) 



where we used that x s — (d — 2 + r/)/2. Specifying T 
as Bk, k = 1, 2, 3, we obtain the (as yet unrenormalized) 
operators 



Bx-H' = -2{Ai+A 2 )-3gA 3 , 
B 2 -H' = -2(A 2 + Ai) - 3gA 5 , 
B 3 -H' = 2(Ai + 2A 2 + Ai) 

+ 6g(A 3 + A 5 ) - 6g 2 (A 6 - A 7 ) . 



(4.11a) 
(4.11b) 

(4.11c) 



These combinations are eigenoperators of the RG, at 
least at O-loop order. At higher loop orders, the renor- 
malization might modify the combinations appearing on 
the right hand sides of Eqs. I|4.11|l . Our explicit calcula- 
tions presented in Sees. IVI and I VII reveals however, that 
this effect is absent at 1-loop order. More importantly, 
general RG arguments (cf. Ref . [24| ) guarantee that this 
effect has no impact on the scaling dimensions at any 
loop order. From Eq. I|4.1U|I we readily deduce that 



XBx-H' 

XB 2 -H' 
XB 3 -H> 



v 



d 

d- 



(4.12a) 

(4.12b) 
(4.12c) 
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Two points are worth being emphasized: (i) the scaling 
relations (|4.12ll are correct to arbitrary order in pertur- 
bation theory, (ii) Equations (|4.11(1 and (|4.12(1 allow us 
consistency checks of our explicit diagrammatic calcula- 
tions. This is particularly valuable because these calcula- 
tions involve fairly many diagrams and one is confronted 
with a certain risk of algebraic errors or erroneous sym- 
metry factors. 

To prepare the ground for the announced consistency 
checks we now evaluate Eqs. (|4.12[) to 1-loop order. Ex- 
ploiting the known e expansion results for </> and v we find 
that the scaling dimensions of the eigenoperators B\ ■ Ji! 
and B3 ■ H.' are given to the order we are working here by 



4e 

XBt-w =d + 2- — + 0(e 2 ). 
x Bi w = d + 2+^- + 0(e 2 ). 



(4.13a) 
(4.13b) 



It is a nice little exercise to derive the result l|4.13bj) ex- 
plicitly in a 1-loop calculation. An insertion of B3 into 
the two-point vertex function I^ 2 ) leads to 1-loop order 
to the singular part 



f ( ^ 2) ({p = o}) = -i 



2/t 



(4.14) 



of the bare function. Thus, it arises a factor 2 in the 
singular part in contrast to a corresponding insertion of 
Bq. Renormalizing B3 = Zg 3 B3 we find, using Z = 1+gj, 
that Zg 3 = 1 — iJ^. From this result we eventually get 



lie 
~21 



m 3 = — + o(s 2 

which then leads to Eq. 14.13b|l . 



(4.15) 



V. RENORMALIZATION OF SPECIFIC 
OPERATORS 

A. RRN specific operators 

To determine the renormalizations of the operators 
Ao , . . . , A3 we study their insertions into the vertex func- 
tions r' n '. We proceed in the spirit of our previous work 
on the RRN. We use 

G( P ,A) = (^(p)^3;(-p))( trunc ) (5.1) 

as the Gaussian propagator, where ip^ip) stands for the 
Fourier transform defined via 



V x (x)= / e^ x Vx(p)- 



(5.2) 



Here, / is an abbreviation for (2tt) d / 2 f d d p. 

(■ • •)( trurLC ) denotes averaging with respect to the Gaus- 
sian part of TL and all fi = 0. Explicitly, our propagator 



reads 



G(p,A) 



5 A,0 



P 2 + wX 2 + T 



jX 2 



(5.3) 

Note that on the right hand side of Eq. (|5.3|1 the propaga- 
tor is decomposed into two parts: one with unrestricted 
values of A (conducting) and one with A = (insulat- 
ing). The decomposition of the propagator leads to a 
decomposition of each of the original (bold) Feynman 
diagrams into an assembly of diagrams made of con- 
ducting and insulating propagators. These conducting 
diagrams resemble essential features of real resistor net- 
works. One can say that they have a real-world interpre- 
tation [E S S, 0, H3, EE EEEE 113, EH 1^, (2^. 

Upon recasting the conducting and insulating propaga- 
tors in the Schwinger representation and using the con- 
tinuum limit in the replica space with D — ► 0, it is easy 
to calculate all required diagrams, see Fig. ^ m dimen- 
sional regularization, at least to one-loop order. For con- 
venience, we work with the rescaled operators 



A' = 



fei/2 



(5.4) 



instead of the original Ai . The fc, are ko — k\ = hi = 
and &3 = 1. The benefit of this rescaling is that the 
primed operators all have the same naive dimension, viz. 
A\ ~ n d+2 . Our 1-loop calculation sketched in Ap- 
pendix [CI gives for the singular parts of the primitively 

(n) 

divergent vertex functions with operator insertions T_^, , 
n = 2 and 3, the results 



r^( P i,p 2 ,A) 



r^(pi, P2 ,A) = wx 2 



wX 2 p 2 [ 1 



5u 

g7 



wX 2 1 



u f P 2 , Pi+pf 
30 



-r^(pi,p 2 ,A) 



-r^( Pl , P2 ,A) 



e \90 
wX 2 



(5.5a) 



(5.5b) 



r Pi+Pi f 1 3u 



2 V lOe 

u (p 2 19wX 2 
1 V 90 + 10 



(5.5c) 



-wX z 



u 1 ' 2 /52wX 2 
45 



.P1+P2 



15 



where p = Pi+p 2 , and 

-r^({ P },{X}) = o, 
-r^({ P },{A» = -- 
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FIG. 1: Diagrams contributing at 1-loop order to the renor- 
malization of the RRN specific operators. These diagrams are 
obtained by first decomposing the bold 1-loop diagrams into 
their conducting diagrams and by then inserting the irrelevant 
operators in the appropriate places. The bold lines symbol- 
ize bold propagators, the light lines stand for conducting and 
the dashed lines for insulating propagators. The labels A and 
B refer to conducting diagrams before insertions are made 
and will be used in the appendix. The dots with wiggly lines 
stand for an insertion of (0) A' , (1) A[, (2) A' 2 and (3) A! 3 , 
respectively. 



In writing Eqs. (|5.5|l and l|5.6l) we have dropped incon- 
sequential factors (r//x 2 )" e/2 and (r/^ 2 )" e / 2 (^ £ /G £ ) 1/2 , 
respectively, for notational simplicity. 

For the purpose of renormalization it is handy to col- 
lect the operators in a vector 



(5.7) 



Due to the mixing, a proper renormalization requires an 
entire renormalization matrix Z. We set 



A'^A 



A' = Z A . 



(5.8) 



This means in turn that the vertex functions with an 
insertion of A!_ are renormalized by 

r (n) . j^(ri) r (n) _ 7 n/2 7 -p(n) 



(5.9) 



The e-poles are eliminated by minimal subtraction. To 
this end, we introduce the Laurent expansion of the Z_- 
matrix 



(5.10) 



Our 1-loop calculation leads to 
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(5.11) 



Bare, unrenormalized quantities are independent of the 
external length scale parameter /i. Thus, we have the 
identity 



M^"^ ({x,\};t,u,w,(j) =0 



(5.12) 



which translates via the renormalizations i|3.1|) and l|5.9|) 
into the RG equation 

d d d d n 

xT%) ({x,A};r,u, W ,/i) =0. (5.13) 
The matrix 7 is given by 
d 



7 



1 fc + u |_ M (i) + if feM (i) 
2= du= 2 V== 



M (1) k 



(5.14) 



Here, k is the diagonal matrix with the diagonal elements 
&0i ki, tt2, and k 3 . To obtain a fixed point solution to 
the RG equation (|5.13|) we recast 7* =7(14*) in terms of 

its eigenvalues and eigenvectors, 



3 

1* =Y1 \ m ) 7 lm(m\ , 

m=0 



(5.15) 



where r\ m are the eigenvalues and (m\ and |m) are the 
corresponding left and right eigenvectors. An RG equa- 
tion for an eigenoperator 



(m\A' 



(5.16) 



is then readily derived form the so-obtained fixed point 
version of Eq. 15.13fl by multiplication with (m|, 



d_ 



d d n 

K T— 1- C W- — T] + 1]„ 



dr 



dw 



xT$l ({ X ,X};r,u*,w,t^ =0. (5.17) 

Using the method of characteristics, it is straightforward 
to solve this RG equation. Augmenting its solution with 
a dimensional analysis to account for naive dimensions 
we find the scaling form 



F 



/2+x„ 



x r%l ({ix, \\ ■ r 1/u t, u\r^ tv w , ^ (5.18) 



with the scaling dimensions x m of the eigenoperator O m 
given by 

x m = d + 2 + r) m . (5.19) 

At this stage our scaling solution (|5.18(l is rather for- 
mal, and we still have to determine the eigenvalues of 7* . 
From Eqs. (|5.11|l and (|B. 14|) we obtain 
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(5.20) 
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At the infrared stable fixed point u* this matrix has the 
eigenvalues 



m 

?/3 



0, 



C* 2 . 

v 



23 



WHS) ^ + otf), 

(23 + 8V30) I ^ + O( £ 2 ). 



(5.21a) 
(5.21b) 

(5.21c) 
(5.21d) 



The left eigenvector belonging to rjo reads (0| = 
(1,0,0,0), i.e., Ao is an eigenoperator of the RG and 
its scaling behavior is governed by r] . From Sec. |W] 
we know rigorously that the scaling dimension of Ao is 
xo = d + 2. Hence, Eq. I|5.21a|) holds to arbitrary order 
in e expansion, r/i is associated with the left eigenvec- 
tor (1| = (0,2,2,3(u*) 1 / 2 ). Thus, we can identify its 
eigenoperator with Bi-H', cf. Sec. lIVI and conclude that 
Eq. (|5.21b|) is valid to arbitrary order in s expansion. To 
1-loop order rji is given by rji = — |f + 0(e 2 ). 

We observe that our 1-loop results for t/q and rji are 
in full agreement with our non-perturbative results de- 
duced from the equation of motion. In other words, our 
1-loop calculation fulfills important stringent consistency 
checks. Of course, these checks do not guarantee the cor- 
rectness of our results for 772 and 773. However, they reas- 
sure us that important features of our calculations like, 
e.g., symmetry factors, are correct. In this indirect sense, 
the consistency checks are in favor of our results for 772 
and 773. 

To extract the sought-after corrections to scaling we 
need to know the scaling behavior of the coupling con- 
stants associated with the RG eigenoperators. The renor- 
malized contribution of the eigenoperators to the Hamil- 
tonian H reads J d d x ^ m v m O m . The flow of the cou- 
pling constants v m under renormalization is therefore de- 
scribed by v m (l) — l u)m v m with the correction to scaling 
exponents 



Vn 



(5.22) 



After all, we find these exponents stemming from the 
RRN specific irrelevant operators to be given by 



u 3 



= <P/v, 

= 2-0.198e + O(e 2 ). 
= 2 + 0.636e + O(e 2 ). 



(5.23a) 
(5.23b) 
(5.23c) 
(5.23d) 



Before we turn to the general percolation operators, 
we find it instructive to briefly re-analyze the work on 
corrections to scaling in RRN by HL. If one erroneously 
neglects the coupling of the operator A\ to the other 
operators under renormalization one is led to a "scalar" 
renormalization factor Zhl that corresponds to the ma- 
trix element Z\ \ of Z_. Thus, one would find to 1-loop 



order 
Ai 
Zhl 

1HL 



Ai — Z H l L A\ , 



-KI- 



0(u 2 



u 
2e 



0(u 2 



~ + 0(u 2 ). 



(5.24) 
(5.25) 
(5.26) 



The anomalous dimension were then t]hl = Ihl* = 
— e/7 + 0(e 2 ), which results in an exponent luhl = 
2 + rjHL = 2 — e/7 + 0(e 2 ). The corresponding crossover 
exponent <p2 defined by HL follows as 02 = 20 — vujhl = 
(2— |f )~ 1 [2(2— |f ) — (2 — f )] = l + 0(e 2 ). Our analysis 
above shows clearly that this exponent has no meaning 
as far as corrections to scaling in RRN are concerned. 



B. Spin specific operators 

To analyze corrections to scaling in continuous spin 
systems, we have to take into account the operator A c 
defined in Eq. (|2.14|l . This operator breaks the 0{D) 
symmetry in replica space. In this sense, its symme- 
try is lower than that of the resistor specific operators. 
This lower symmetry has an important consequence. 
Though the renormalization of A c generates the oper- 
ators ^Iq, . . . , A' 3 , these do not in turn generate A c under 
renormalization. 

This structure is very similar to the one we encountered 
in studying multifractality in RRN [ltA |20| and random 
resistor diode networks |2ll |22| . In the field theoretic 
description of these networks, there are dangerous ir- 
relevant operators corresponding to the multifractal mo- 
ments of the current distribution on the networks. These 
operators (masters) generate a whole bunch of other op- 
erators (servants). The servants on the other hand, do 
not generate their masters. All servants must be taken 
into account in the renormalization process, at least in 
principle. However, the renormalization matrices have a 
particular, simple structure. Due to this simple struc- 
ture, the scaling index of a master operator is completely 
determined by a single element of the renormalization 
matrix. Hence, for the practical purpose of calculating a 
masters scaling index, the servants can be neglected. 

To facilitate the renormalization of the spin specific 
operators, we introduce the vector 



Alss — (^c, A , . . . , .A3)" 



(5.27) 



A proper renormalization requires a 5 x 5 renormalization 
matrix Z_ which we introduce by setting 



Aw * A* 



A a = Z A a 



(5.28) 



The arguments given above imply that Z_ is of the form 



Z c Xc 



(5.29) 



where Z_ is the renormalization matrix defined in 
Eq. (|EH^ and = (0,0,0,0). The elements of 2L C must 
be chosen, in principle, to cancel e poles associated with 
the servants generated by A c . In practice, however, we do 
not need to determine these elements for computing the 
correction to scaling exponents. In a 1-loop calculation 
we find 



(4) 



Z r = l 



17 u 9 . 



(5.30) 



The matrix 7 , which is given up to an obvious mod- 

ification by Eq. I]5.14[l . inherits the simple structure of 
Z_ . As a consequence, the eigenvalues of 7 at the fixed 

point u* are 770 to 773 and rj c . From Eq. (|5.30|1 we obtain 
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This eigenvalue leads finally to a correction exponent 



(5.32) 



in addition to luq- 
HL is related to u> c via 
1-loop order 



, UJ3 . The crossover exponent ip c of 
= 2<fi — voj c . We obtain to 
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T55 + < e2) 



(5.33) 



in conformity with the value given by HL. 
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FIG. 2: Diagrams contributing at 1-loop order to the renor- 
malization of the general percolation operators. The bold 
lines symbolize bold propagators. Single dots accompanied 
by a wiggly line respectively stand for an insertion of (4) A4, 
(5) A'5, and (7) Aj. An insertion of A'q is depicted with help 
of 2 dots and 2 wiggly lines. Note that we have, in contrast to 
Fig- displayed bold instead of conducting diagrams to save 
space. 



VI. RENORMALIZATION OF GENERAL 
PERCOLATION OPERATORS 

Now we turn to the renormalization of the operators 
A4 to At- The corrections to scaling arising from these 
operators have been calculated a long time ago by Amit et 
al. [2^|. Nevertheless, we think that the general percola- 
tion operators deserve some attention here for 2 reasons: 
(i) Amit et al. used the usual Potts model formulation 
of percolation that is different from the RRN formula- 
tion in the way constraints on the order parameter field 

are implemented, viz. X)g s ( x i^) = ^ versus V'o( x ) = 0- 
It seems desirable to have a treatment of both the spe- 
cific and the general percolation operators that is self- 
contained within one formulation, (ii) Because the calcu- 
lation is somewhat involved, two independent approaches 
help to guarantee the correctness of results. In fact, our 
results concerning the general percolatio n op erators co- 
incide in the end with the results of Ref . 28] . 

The general percolation operators affect the renormal- 
ization of the 2, 3 and 4-point vertex functions. To asses 
this effect, we have to calculate these vertex functions 
with insertions of A^ to Aj. The Feynman diagrams 
that contribute to these vertex functions are summarized 
in Fig. |21 To calculate these diagrams, we decompose 
them into their conducting diagrams and then take the 
limit w —* 0. Then, the conducting and the insulating 



propagators can be replaced by simple l/(r + p 2 ) prop- 
agators. For most of the diagrams, decomposition cul- 
minates into to a simple numerical factor. A little cau- 
tion must be exercised, however, to distinguish between 
terms corresponding to Aq and A7, respectively. To fos- 
ter this distinction, we keep the dependence of these two 
operators on the replica currents explicit. This depen- 
dence is different in the way the sum over the external 
replica currents, say Ai, • • • , A4 vanishes, viz. Ae is pro- 
portional to S = {S Xl+ x 2 ,o a A 3 +A 4 ,o + S x 1 +x 3 M S x 2 +x i M + 
^a"i+a 4 o^a 2 +a 3 0)/^ wnereas A7 is proportional to F = 

^Ai + A 2 +A 3 + A 4 ,0- 

As for the specific operators, it is convenient to rescale 
the general operators. To be specific, we rescale the gen- 
eral percolation operators according to Eq. I|5.4|l with 
ki = 0, fcs = 1 and k$ = kj = 2. 

Our 1-loop calculation gives for the 2-point functions 
with insertions 
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where we have dropped inconsequential factors 
(r//x 2 )~ e / 2 for notational simplicity. For the 3-point 
functions with insertions we find 
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Here, we omitted factors (t//j 2 ) £ / 2 (^ e /G e ) 1 / 2 . Our re- 
sults for the 4-point functions with insertions read 
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(6.3a) 
(6.3b) 
(6.3c) 
(6.3d) 



where we dropped factors (r / ' p?)~ e l 2 'fj, e / 'G e . 

The remaining RG analysis proceeds completely anal- 
ogous to that for the RRN specific operators. The RG 
equation governing the general percolation operators has 
the same structure as Eq. (I5.13|l . Of course, here we can 
set w — 0. Basically, we just need to make the replace- 
ments 
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vAgp — (^4 j ■ • • ! <A<j) , 



(6.4) 
(6.5) 



with the latter matrix resulting from 1)6. lfl . I|6.2|) and (|6.3|) 
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At the infrared stable fixed point, 7 has the eigenvalues 

=gp 
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= 26 



/889)^- + 0(£ 2 ). 



(6.7a) 
(6.7b) 

(6.7c) 
(6.7d) 



The eigenvalues 774 and 775 are associated with the 
left eigenvectors (4| = (2, 3(u*) 1 / 2 , 0, 0) and (5| = 
(1, 3(m*) 1 / 2 , — 3u*, 3u*), respectively. By comparison 
with the w — > limit of Eqs. (|4.11l) . we identify B 2 ■ H' 
and 63 • H' as the corresponding eigenoperators. It fol- 
lows that Eqs. I|6.7al) and 1)6. 7b|) are correct to arbitrary 
order in perturbation theory. The eigenvalues stated in 
Eqs. (|6.7|) entail the correction to scaling exponents 



U4, = 2, 

uj 5 = 2 — 77 , 

lu 6 = 2-0.182e + O(e 2 ), 

lu 7 = 2 + 2.658e + 0(e 2 ) 

in full agreement with the results of Amit et 



(6.8a) 
(6.8b) 
(6.8c) 
(6.8d) 



VII. CRITICAL BEHAVIOR OF THE AVERAGE 
RESISTANCE ETC. 

Now we can harvest the results of our RG analysis to 
write down scaling expressions, including the most im- 
portant corrections, for key observables for the RRN and 
the ir-y-model. We will elaborate on the RRN in some 
detail and assemble step by step the generating function 
G(x, x',A) from which we then extract and E. To 

cover spin models, the generating function requires some 
modifications that will be explained. 



A. Random resistor networks 

The generating function has two kinds of ingredients, 
viz. the 2-point correlation function without insertions 
and the set of 2-point correlation functions obtained by 
inserting each of the RG eigenoperators. From Eq. (|3.3|l . 
Eq. (|5.18l) as well as its counterpart for the general per- 
colation operators we collect that 



g(x,x',a) = f*- 2 +"i G< 2 > (^|x-x'|,r^A 2 ) + XS^*Gg ^|x-x'U-^A 2 ) +••• I. (7.1) 
^ fe=i ' 

Here, we used that the 2-point correlation function is the inverse of the 2-point vertex function and we applied Fourier 
transformation to switch from momentum to position space. The ellipsis stands for contribution from irrelevant 
operators with a naive dimension higher than 8. Next we expand Eq. I|7.1|l in a power series in wX 2 as well as in the 
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deviation u — u* from the fixed point. Choosing £ = |x — x'| 1 and setting all non-universal constants equal to 1 for 
notational simplicity we get the leading terms 



g(x,x',AJ = |x- x 'r (d - 2+?,) 
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(7.2) 



By virtue of Eq. I|2.9|) we know that we now can extract 
the average resistance by taking the derivative with re- 
spect to — A 2 /2 and then setting A = 0. After all, we 
arrive at 
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Here, A and A/, are non-universal constants. We opted 
to write them down explicitly here to make closer contact 
to Eq. (|1.4|) . A goes to zero for u — ► u* . The Ak vanish 
for Vk — > 0. 

Next we turn to the average conductivity of the RRN. 
Commonly, the conductivity £ of percolating systems is 
defined with respect to a bus bar geometry where the 
network is placed between two parallel superconducting 
plates (the electrodes) of area L d ~ x a distance L apart. 
From the above we expect that the average conductance 
cr of this finite-size system scales as 



<t(L,t) 



\T\ 
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fe=l 



(7.4) 



where £ ~ |r| "is the percolation correlation length and 
II, 5 and the 2^ are scaling function with the properties 



U{x) ~ E(x) ~ E k (x) 



const for x <C 1 

x d - 2 for i » 1 ' (7 ' 5j 



Above the percolation threshold (t < 0) the RRN be- 
haves on length scales large compared to the correlation 
length £ ~ \t\~ v like a homogeneous system with con- 
ductivity S. Hence, we may write for L>( that 

S(r) ~ L 2 ~ d a(L, t) . (7.6) 

Merging Eqs. f7~l|) . lf73|l and 1(73)1 we finally get 



1 + 5M 



fc=i 



(7.7) 



with non-universal amplitudes B and Of course, B 
vanishes for u — > u* and the bk vanish for Vk — > 0. 



B. Spin systems 

For spin systems like the x-y-model we had to include 
the irrelevant operator A c into our RG analysis because 
these systems are, unlike the RRN, not 0(Z?)-invariant 
in replica space. Thus, the RG has 9 irrelevant eigen- 
operators of naive dimension 8 and the counterpart of 
Eq. H7.1fl has an extra term that features the exponent 
u> c . By applying basically the same steps as for the av- 
erage resistance we find that the first cumulant of the 
angular fluctuations scales as 
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+ C c \x- x'\-"° + °k \x ~ x'\^ k + 



k=l 



where C, C c and the Ck are non-universal amplitudes. 

Equation l|7.8(l concludes our results on corrections 
to scaling. However, our analysis also sheds light on 
the second cumulant of the angular fluctuations. Equa- 
tion (|2.11() implies that we can derive C^p via taking 
the derivative with respect to ^(A). This homogenous 
polynomial in A exclusively appears in A c . We saw that 
A c has the property of being a master operator whose 
scaling behavior is governed by lu c . Consequentially, we 
obtain that 



(7.9) 



Here, A c gives rise to the leading scaling behavior, i.e., 
A c is a dangerous irrelevant operator as far as C(p is 
concerned. 

Of course it is interesting in this context to ask, what 
the leading scaling behavior of the higher cumulants 
might be. It can be shown to arbitrary order in per- 
turbation theory [3(|, that 



Cf ~ Ix-x'l^/" + ■■■ 



(7.10) 



where the ipi are the critical exponents of the multifractal 
moments 



«f - (e (f 
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of the current distribution on RRNs, with if, being the 
current flowing through bond b and I being the to- 
tal external current. Note that the HL result (j> r = 
1 + e/105 + 0(e 2 ) and the e expansion result pLUal^ 
for tl>2 (the entire family of the ipi is known to two- loop 
order [H 20]) are in full agreement. The upshot here is 
that Eq. (|7.9|l represents merely an instance of the gen- 
eral result l|7.10|l . 



VIII. CONCLUDING REMARKS 

We have studied corrections to scaling in RRN and 
continuous spin models. As far as the leading scaling 
behavior and its leading corrections are concerned, the 
HL model provides a unified description of both systems. 
Being interested in next to leading corrections, however, 
we had to consider distinct sets of irrelevant operators 
for the two systems. 

In both systems, we found the typical mixing of irrel- 
evant operators under renormalization. Thus, we had to 
compute an entire renormalization matrices. This is the 
reason why we restricted ourself to considering irrelevant 
operators with a naive dimension 8. At least in principle, 
one could analyze higher corrections to scaling originat- 
ing from irrelevant operators of naive dimension 10, 12 
and so on. However, in these cases the renormalization 
matrices become prohibitively big and their computation 
and diagonalization requires enormous effort. 

One of the spin specific operators, namely A c , is quali- 
tatively different from the remaining irrelevant operators 
under consideration. A c has the properties of a master 
operator. Hence, it is sufficient to calculate a single renor- 
malization constant to determine its scaling dimension. 
This can be done in an el egan t way and with moderate 
effort up to two- loop order [3fJ . One has to bear in mind, 
though, that the corrections to scaling are not only re- 
sulting from A c but also from the other spin specific as 
well as from the general percolation operators. Over all, 
it would require a lot of work to extend our results on cor- 
rections to scaling to 2-loop order or to include irrelevant 
operators of higher naive dimension. 

It is interesting to compare the corrections to scaling 
stemming from the specific and the general percolation 
operators with the corrections arising from the presence 
of a surface. In Ref. Q we studied a semi-infinite RRN 
at the so-called special and ordinary transitions |2|j . We 
calculated the corrections to induced by the sur- 

face when the terminal points x and x' are located on 
the surface. At the special transition we found a cor- 
rection that vanishes as |x — x'|~ ws for increasing ter- 
minal separation with the correction to scaling exponent 
ujg = 1 — e/21 + 0(e 2 ). This correction falls of slower 
than the corrections induced by the specific and the gen- 
eral percolation operators, i.e., this correction represents 
the next to leading correction if a surface is present. At 
the ordinary transition we found a correction |x — x'|~ Ws 
with a/g? = 3 — 23e/105 -I- 0(e 2 ), i.e., this correction van- 



ishes faster than any of the corrections induced by the 
specific and the general percolation operators. 
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APPENDIX A: USEFUL IDENTITIES 

Equation 14.811 can be regarded as the centerpiece of 
Sec. lIVI In this appendix we first derive a Ward identity 
that relates the vertex functions with a single insertion 
of the operator 7i to the vertex functions without inser- 
tion. Then we obtain Eq. I|4.8|l as a corollary by slightly 
modifying our previous arguments. 



1. A Ward identity 

Now we augment the HL model with an external field, 
i.e., we consider the Hamiltonian 

n h =n-(h,s), (Ai) 

with H given by Eq. 1)2. 4|) with /; = and where we have 
used the abbreviated notation 

(h,s) = J d d x J^h(x,6)s{x,6) . (A2) 

To facilitate our argument, we consider the shift 

s(x,9) s(x,0) + c(x,9) (A3) 

of the order parameter field. This shift leads to the mod- 
ification 

H h ^H h + (c,H' h ) (A4) 

of the Hamiltonian. Here and in the following we omit in- 
consequential terms of order 0(c 2 ). By virtue of Eq. (|A4|) 
an application of the shift to 

((h, s) n ) = J Vs (h, s) n exp (-74) (A5) 

yields 

((h,s) n (c,H' h ))=n(h,c) (iKs)* 1 - 1 ) . (A6) 

Dividing both sides of Eq. (|A6() by n\ and summing from 
1 to oo leads to 

Z ic , ni) [h] = (h,c)Z[h], (A7) 
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where 



Z[h] = (ew(h,3))=f;&yp- (A8) 

is the generating functional for the correlation functions 
of the order parameter field and 

f M(c,?4) (M)") rAq , 

n=l 

is the generating functional for the corresponding corre- 
lation functions with an insertion of (c,7i' h ). Defining 



Z[h,c} = Z (c ^Jh} + Z [ h] 



(A10) 



and switching to the generating functional of connected 
correlation functions, 



W[h,c] =]n(Z[h,c]) 
and so on, we find 

W[h,c]=W[h] + (h,c). 



(All) 



(A12) 



Moving over to the generating functional T[h, c] of vertex 
functions via the Legendre transformation 



with 



T[s,c\+W[h,c] = {h,s). 



5T _ 5W _ 5T _ SW 
Ss Sh 5c 5c 



we arrive at 



T[s,c] = (h,s)-W[h] - (h,c) 



(A13) 



(A14) 



(A15) 



Now we take the functional derivative with respect to c 
to obtain the identity 



T w 



ST _ 5T 
5c 5s 



(A16) 



between the generating functional of vertex functions 
with an insertion of TC and the generating functional T of 
the usual vertex functions with out any insertion. Finally, 
we take n functional derivatives with respect to the order 
parameter field to obtain the Ward identity 



r(") _ r(™+ 1 ) 
H' 1 



2. Derivation of Eq. 14, 8[) 



(A17) 



In order to derive Eq. Q4.8p . we just need to modify our 
arguments of Appendix IjA 1(1 slightly. Here, we assume 
that c is a composite field comprising the order parame- 
ter field and possibly its gradients in real and in replica 
space. Consequentially, we have to replace Eq. (|A6(I by 



where we once more have omitted terms of order 0(c 2 ). 
Upon taking n functional derivatives with respect to the 
external field we obtain 

(«(x 1 )...«(x n )(c,W'))= (A19) 

71 

= J~] (a(xi) ■ • ■ s(xj-i) c(xj) s(xj + i) • ■ ■ s(x n )) . 

i=l 

Setting c(x, 0) = a(x, 9)J : (x, 9), where a is independent 
of the order parameter field, and by taking the functional 
derivative with respect to a we obtain 



s(x 1 ,0i)--.«(x n ,0 n ).F(x, 0)W'(x,< 

n 
i=l 

x s(x i+ i, e i+1 ) ■ • • s(x n , 4)) • (A20) 

Equation l|4.8|l is now readily obtained by integrating 
over 9. 



APPENDIX B: COMPOSITE FIELDS IN THE HL- 
AND THE POTTS MODEL 

In Sec. IIVI we exploited information on lower dimen- 
sional operators to draw conclusions on the scaling di- 
mensions of several eight-dimensional operators that are 
associated with corrections to scaling. This appendix is 
intended to provide some background on the RG behav- 
ior of the lower dimensional operators and to establish 
some of the results that serve as an input in Sec. IIVI 

The following arguments refer primarily to the HL 
model. Since the HL model reduces to the usual Potts 
model with N — (2M) D states upon setting w = 0, how- 
ever, our reasonings also apply to the latter model. Being 
interested in the Potts model, one basically just has to 
set w = in any of the formulas in this section. To make 
closer contact to the conventional notation for the Potts 
model, one may replace s(x, 9) by Si(x) with i = 1, • • • , N 
(along with h(x, 9) — > /ij(x) for the external field and so 
on) and the integral fx by a summation over i. 

1. Four-dimensional operators 

We start by considering composite operators with the 
naive dimension four. The simplest operators of this kind 
are Bo, B\ and B^. The scaling dimensions of these op- 
erators follow directly from known RG results. We will 
revisit these operators briefly towards the end of this sub- 
section. 

Bo belongs to the trivial representation of the permu- 
tation symmetry group. Its counterpart belonging to the 
fundamental representation is 



((h,s) n (c,H' h ))=n((h,c) (M) B_1 ) 



(A18) 



Cl = s ~N 



(Bl) 
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The scaling dimension of this operator does not follow 
immediately from known results. It will be derived in 
the following. We consider the Hamiltonian 

n a = |^^|^ + I^ + i (Vs)2 + | (Ves)2 

(B2) 



I- 3 



hs 



where h = /i(x, 9) is an external field and a = c(x, 9) 
has the property Jga = 0. Due to the extra terms, we 
need renormalizations in addition to those specified in 
the renormalization scheme (|3 -II) . viz., 



a = Z Z„a . 



(B3a) 



h -> h = Z~ 1/2 (h + Act 2 + Bra + CV 2 a + DwW 2 g aj . 

(B3b) 

In order to determine the renormalization factor Z a and 
the additive renormalization constants A to D, we per- 
form a shift 

s = s' +c, (B4) 

where c is assumed to satisfy the condition Jgc — 0. This 
shift transforms the Hamiltonian so that 



Ha [s, t, a, w, g, h] = H (T [s',T,a',w,g, h'] 
+ H a [c, t, a,w,g,h] , 



where 



a = a + gc , 

h! = h — tc — ac + V 2 c + wV«c - 



(B5) 

(B6a) 
(B6b) 



In other words: the Hamiltonian TL a is form invariant 
under the shift i|B4fl up to an inconsequential term that 
does not depend on the order parameter field. Note that 
the shift neither modifies r nor w or g. Since the Z 
factors depend only on the dimensionless variant u of g, 
it follows that none of the Z factors is affected by the 
shift. 

Equation (|B4fl implies that the renormalized version of 
c is given by 

c = Z- 1/2 c. (B7) 
Exploiting this we find by renormalizing Eq. ljB6afl that 
Z a = Z g . (B8) 



ti = h- (2gA + Z g )ac - (2gA + Z g )^c 2 - (Bg + Z t )tc 



Renormalization of Eq. (|B6b|) yields 



(Cg + Z)V 2 c - (Dg - Z w )wV 2 g c . 



(B9) 



Demanding that the renormalized hi retains its original 
from, cf. Eq. l|B6bJl . we obtain 



(BlOa) 
(BlOb) 



9 

Z w - 1 



25 ' 
9 



9 



for the additive renormalization constants. Now we know 
all the renormalization factors and constants and hence 
we can write down the renormalized version of the Hamil- 
tonian Ti a , viz. 



where 



Ha =H + {a, A) +0(a 2 ) , 



(Bll) 



TS 



Z - 1 



V 2 s 



Z,„ - 1 



9 



9 



(B12) 



is the fully renormalized version of the operator T>i. 
Naively, one might have expected that a mere multiplica- 
tion with Z^ 1 were sufficient to renormalize T>±. Equa- 
tion l|B12|l . however, shows clearly that this is not the 
case. 

Next we determine the scaling dimension of A via 
the scaling behavior of its coupling a. Upon taking the 
derivative of Eq. (|B3a|) with respect to the external in- 
verse length scale fi we obtain 

k = fid^alo = 7 - lg ■ ( B1 3) 
Recalling that 7* = 77 and deducing from 

= (-£ + 3 7 - 2 7g )u (B14) 
that 7* = — e/2 + 3ry/2 we find 

«* = ^ (B15) 

as the fixed point value of k. Taking into account that 
the naive dimension of a is two, we get that the RG flow 
of a in the vicinity of the fixed point is given by 



with 



a{t) =or v ° 
d - 2 + rj 



Finally, this leads to 



xa = d- y„ = 



d + 2-r) 



(B16) 



(B17) 



(B18) 



for the scaling dimension of A. Note that x_a is identical 
to the scaling dimension of the operator B3. 
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We have to point out that A is not an eigenoperator 
of the RG. This fact follows from two observations: (i) 
$3 = TC is an eigenoperator, as can be seen from the 
Ward identity HA17(I . and (ii) A is just one ingredient of 
B 3 , 



#3 = gA + rs — V 2 ,s — wVgS , 



(B19) 



where, of course, all quantities in this equation are renor- 
malized quantities. For completeness we mention that 
the Ward identity (|A17|) takes on the form 



^r!; ( ) p) ({q}) = -r(™+ 1 )({q},p) 



(p 2 + wX 2 



r) 6 n ,i ■ 
(B20) 

when expressed in terms of A. 

As announced above, we now briefly return to the sim- 
pler four-dimensional operators. The scaling dimension 
of Bq follows immediately from the scaling behavior of 
its coupling constant t(£) = t£~ Vt with y T = 1/v, 



x B a = d — y T = d — 1/v . 



(B21) 



The operators B\ and B2 are obtained by applying, re- 
spectively, wV 2 , and V 2 to the fundamental field s and 
hence 



d+2 



x Bl = x s 
x@ 2 x s 



V 



d-2 + r) <j) 
2 + v 



(B22a) 
(B22b) 



By now, we have expressed the scaling dimensions of all 
the four-dimensional operators which enter our analysis 
in Sec. II VI in terms of the fundamental exponents 77, v, 
and d>. 



2. Six-dimensional operators 

The most self-evident operators with naive dimension 
six are, perhaps, 



C2 = wsVgS , C3 = sV s , C4 



(B23) 



One eigen-combination of these operators follows imme- 
diately from Eq. JOJ, viz. s ■ TC. Due to Eq. (|4~TU|) we 
know that the scaling dimension of this eigenoperator is 



Xs-W 



(B24) 



i.e., it is marginal in any dimension. The physical reason 
for this marginality is that s ■ TC can be removed by a 
rescaling of the amplitude of the order parameter field. 

Next, we take a closer look at the scaling dimen- 
sion of the renormalized version of D4. Setting all the 



fi in Eq. (|2.4|l to zero, applying the renormalization 
scheme l|3.1|l and taking the derivative with respect to 
the renormalized coupling constant g we obtain 



where 



-T- = I d a x 

dg 



B. 



(B25) 



B 



\[z a + 2uZ' g y + 



^( Vs ) 2 + ^(v 9S ) 2 
9 9 



uZ' T 
9 



-TS 



(B26) 



is the fully renormalized version V4. Here and in the 
following the prime indicates derivatives with respect to 
u. Equation (|B25I) implies the Ward identity 



dg 



_ -p(n) 



(B27) 



Taking into account that the genuine independent vari- 
ables of the RGE for the vertex functions 1^™) are fi, u, 
t, and w, and exploiting the commutator 



d_ 

dg 



e d 
Ydg 



(B28) 



we find that the RGE for the vertex functions with an 
insertion of B is given by 



(jid^ + kt8 t + (wd-u, + (3d u ~ -^1 + P' + -j 



(n) 



-{nj 
9 



2k't<9 t - 2C'wd w )T 



(B29) 



Note that this RGE is not homogeneous and hence B is 
not an eigenoperator. Nevertheless, we can deduce the 
scaling dimension of B from the homogeneous part of the 
RGE (|B29|) . Taking into account the operator's naive 
dimension 3(d — 2)/2, we obtain 

x B = d + u, (B30) 

where u> is the Wegner exponent featured in Sec. IIIIBI 



APPENDIX C: CALCULATION OF FEYNMAN 
DIAGRAMS 

In this appendix we give some details on the diagram- 
matic calculation that leads to Eqs. 15.5fl and l|5.6|l . In- 
stead of elaborating on all of the diagrams we restrict 
ourself to a couple of representative examples. The steps 
explained at these instances can then easily be adapted 
for the remaining diagrams. 



J 



We start by considering diagram A with an insertion of Ai, cf. Fig.^ Upon writing the conducting propagators in 
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Schwinger representation we have 



A Al = g 2 W 



/>oo 

/ dsids2ds 3 exp [— (si + s 2 + s 3 )t] 
Jo 

I Yl (« 2 ) 2 ex P { - s i I k2 + wk 2 }~ s 2 [(k + p) 2 + wk 2 ] - s 3 [(k - q) 2 + w{k - A) 2 ] } 



(Cl) 



Note that symmetry factor of A - 4 1 is 1 and not 1/2 (the symmetry factor of A), because there are 2 possibilities (the 
2 conducting propagators) to insert A\. Now we carry out a completion of squares for the momenta as well as for the 
currents. After the straightforward momentum integration we arrive at 



2„,,2 



^Ai = 9 w 



x exp 



ds\ds 2 ds 3 



(4tt)<V 2 J Q ( Sl + S2 + S3 )d/2 



-R(s 1 ,s 2 ,s 3 )w\ 2 \ ^2 



exp [-(si + s 2 + s 3 )r] exp 



S2S 3 (p + q) 2 + sis 2 p 2 + sis 3 q 2 



«1 + S 2 + S3 



S3 



Sl + s 2 + s 3 



exp [-(si + s 2 + s 3 )k 2 ] 



(C2) 



where i?(si, s 2 , s 3 ) = (sis 2 + s 2 S3)/(si + s 2 + s 3 ) is, according to our real-world interpretation, the total resistance of 
the diagram A Al . At this stage it is useful to switch to continuous loop currents and to replace the summation 
by the integration By standard Gaussian integration we then find for D — > 



2„,,2 



ds\ds 2 ds 3 



(471)^/2 J Q ( Sl + S 2 + S 3 ) d / 2 



exp [-(si + s 2 + s 3 )r] 



(si + s 2 + s 3 ) 



A 2 



A' 2 



(si + s 2 + s 3 ) 4 w 



[s 2 s 3 (p + q) 2 + sis 2 p 2 + sis 3 q 2 + (sis 3 + s 2 s 3 )wA 2 ] } , 



(si + s 2 + S3) 3 w 

(C3) 



where we have carried out a Taylor expansion of the second and third exponential function appearing in Eq. I|C2|) 
and where we have discarded all convergent terms. The remaining integrations over the Schwinger parameters can be 
simplified by setting si = ta, s 2 = ty and S3 = t(l — x — y) and then integrating t from to 00, y from to 1 — x and 
x from to 1. Expanding the so obtained intermediate result for small e, we obtain 



-g 2 wX- 



Ge 



r -e/2\r P 2 P 2 + P 2 ] 

1 3 90 30 J ' 



(C4) 



where we have set p = p 1 + p 2 and q = — p 2 . The computation of B Al is, in comparison, simple because it does not 



i 2 \ 

involve a summation over a loop current. The total 1-loop contribution to — T^; is given by 



(C5) 



Here we have a clear example of the mixing of irrelevant operators in our perturbation calculation. The insertion of 
A\ does not only generate primitive divergences proportional to w 2 (A 2 ) 2 but also those of the type twX 2 and p 2 wX 2 . 
Note the we dropped the first term in the braces in Eq. I|5.5|l because it leads only to subdominant correction due to 
the factor r. 

As a next example, we consider the diagram A^ 2 that stands for 



k 2 + (k + p) : 



A^ 2 = g 2 w / dsids 2 ds 3 exp [— (si + s 2 + s 3 )r] / Vj it 
JO Jk z 

x exp { - si [k 2 + wk 2 ] - s 2 [(k + p) 2 + wk 2 } - s 3 [(k - q) 2 + w(k - A) 2 ] } 

Completion of squares in the momenta and the currents leads to 



(C6) 



g 2 f°° 

Aa 2 — — w / ds 1 ds 2 ds 3 exp [— (si + s 2 + s 3 )r] exp 
^ Jo 

-i?(si, s 2 , s 3 )w;A 2 



S2S 3 (p + q) 2 + sis 2 p 2 + sis 3 q 2 



x exp 



S3 

si + s 2 + s 3 



Sl + s 2 + S3 

exp [— (si + s 2 + s 3 )k 2 ] 



s 2 p - s 3 q 

Sl + s 2 + S3 



k + p 



s 2 p - s 3 q 

Sl + s 2 + S3 



exp [-(s x + s 2 + s 3 )k 2 ] , 



(C7) 
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where we switched to continuous loop currents. Carrying out both Gaussian integrations and performing a Taylor 
expansion we obtain in the replica limit 



A A2 = ^wX 2 



(4 



i r 

Jo 



dsids2ds3 
(Si + S 2 + s 3 y/ 2 



exp [-(si + s 2 + s 3 )t] 



+ 



*§P 2 



(si + s 2 + s 3 ) 3 (si + s 2 + S3) 2 



2 sjs 2 p 2 - 2s 2 s 3 p ■ q + s 4 q 2 _ 2 s 2 s 2 p 2 - s 3 p ■ q _ rf s 2 s|(p + q) 2 + s lS2 s 2 p 2 + SiS^q 2 



d 



(Sl +3 2 + S3) 4 

(sis| + s 2 sl)wX 2 ' 



(si + s 2 + S3) 3 



(si + s 2 + S3) 4 



(C8) 

(si + s 2 + s 3 ) 4 1 

where we have discarded convergent terms. Using the same change of variables as for A Al we integrate out the 
Schwinger parameters. This leads in e expansion to the result 

ro Ge 

e 



K M =-g 2 wX 2 — T -e/2lL 



p' , p? + pi 



10 90 



60 



(C9) 



(2) 

The diagram B^ 2 can be calculated by similar means. For the entire 1-loop contribution to — P^; we find 

- r ( 2 ) i-ioop = A 2 B = f w y?9l T -/2 ( 5r + + E! + ^Ei±P|) \ . 

-42 .y £ \ 2 10 90 20 J 

Diagram A Aa serves as our final example. In Schwinger representation this diagram reads 

r a a r f ^ 1 fsr « 2 + (^-A) 2 + A 2 

A^ 3 = -510 y dsids 2 exp [-(si + s 2 )rj / 2^ 



(C10) 



k s 



x exp { - si [(k - q) 2 + w(it - A) 2 ] - s 2 [(k + p) 



2 1 -2 

1 10K 



]}• 



(Cll) 



Once more we undertake a completion of squares in the momenta and currents. After integrating out the loop 
momentum we have 



A^ 3 



9 1 
■ w ■ 



dsids2 



3 (4n) d / 2 J (si + s 2 ) d ' 2 

x / [a 2 + U-^^ x]\ U 

JkY v si + s 2 / V 



exp [-(si + s 2 )t] exp -R{s ll s 2 ) [wX 2 + (p + q) 2 ] 



si 



-A 



exp [-(si + s 2 )i 



(C12) 



/ L "-' - Sl -|- S 2 

with the total resistance of this diagram being R(si, s 2 ) = sis 2 /(si + s 2 ). Integration over the loop current and Taylor 
expansion gives, up to convergent terms, 



M 3 



-^wX 
3 (4 



— f 

*) d/2 Jo 



ds\ds2 

{81 + S 2 ) d ' 2 



exp [-(si + s 2 )r] I I + 2 



(si + s 2 ) 2 



(C13) 



- R( Sl ,s 2 ) [wX 2 + (p + q) 2 ] - 2 - Y 2 [ W A 2 + (p + q) 2 ] 

(si+s 2 ) J 

The integrations over the Schwinger parameters can here to simplified by setting si = tx and s 2 = t(l — x). After 
symmetrizing the external momenta the e expanded result reads 



A^ 3 = g wX —t 



-e/2 



2 fl0r 8wA 2 8(p 2 + p^) 



9 + 45 



90 



The calculation of B^ 3 is fairly easy. Merging the 2 results we obtain 



(2) l-loop 



Ge 



A M -2B As = - gw X 2 —T-^ 2 \ — 



14 t 52 wX 2 2(p 2 +p 2 .) 



45 



15 



(C14) 



(C15) 



The 3-leg diagram with insertions can be computed by similar techniques as we have used for the 2-leg diagrams. 
The 3-leg diagram benefit from the extra simplification that they can be evaluated at vanishing external momenta 
because all their contributions with non zero external mome nta are convergent. 
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